Spin Precession and Avalanches 
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In many magnetic materials, spin dynamics at short times are dominated by precessional motion 
as damping is relatively small. In the limit of no damping and no thermal noise, we show that 
for a large enough initial instability, an avalanche can transition to an ergodic phase where the 
state is equivalent to one at finite temperature, often above that for ferromagnetic ordering. This 
dynamical nucleation phenomenon is analyzed theoretically. For small finite damping the high 
temperature growth front becomes spread out over a large region. The implications for real materials 
are discussed. 
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Ferromagnetic systems that are subject to slowly 
changing external magnetic fields very commonly show 
avalanche-like responses i^. This leads to hysteresis, as 
avalanches occur over a very fast timescale resulting in 
irreversibility and entropy production. 

A large amount of experimental and theoretical work 
has been devoted to understanding aspects of this be- 
havior, such as Barkhausen noise [l| which demonstrates 
that there is often a large degree of reproducibility in 
the mesoscopic dynamics on repeated cycling of the field 
and also interesting critical properties 0, 1^ ■ With ad- 
vances in experimental techniques, direct tests of the re- 
producibility of magnetic memory have been undertaken 
recently which highlighted the prominent role of sample 
disorder |j]. 

The theoretical treatments to date have largely relied 
on simplified models such as the Ising model to under- 
stand these complex systems. The dynamics of such 
models have been purely relaxational, the extreme limit 
of large damping, whereby a spin is flipped if the energy 
of the system is decreased by doing so, the excess energy 
being transferred out of spin degrees of freedom, to for 
example, phonons. These models have had a great deal 
of success in describing many features of disorder fer- 
romagnets, showing fascinating properties, for instance 
"Return Point Memory" (RPM) dQ. 

However real magnets are typically dominated by pre- 
cessional effects on short enough time scales. The 
Landau-Lifshitz- Gilbert (LLG) equation [3, contains a 
precessional term and a dissipative one 



ds ,_^ ^. 

— - = -s X (B-7S X B), 
at 



(1) 



where s is a microscopic magnetic moment, B is the lo- 
cal effective field, and 7 is a damping coefficient. 7 mea- 
sures the relative importance of damping to precession. It 
ranges from about .01 to 1 in real materials 6J. Therefore 
in many magnetic materials, there should be an interest- 
ing short time regime where it makes sense to regard the 
damping as a perturbation. In fact, we will show that 
there is a short time scale in which the magnetic response 
to the applied field change is strongly influenced by the 



finite level of damping in real materials. Furthermore, 
even macroscopic and long-time properties, such as hys- 
teresis loops are influenced, by the level of damping and 
other materials properties, that we shall explain in more 
detail below. 

Throughout this paper, we only consider the case of 
zero thermal noise. This is because we will see that ef- 
fective finite temperature behavior is found even in this 
case and we want to carefully separate out these two ef- 
fects. 

If we adiabatically lower the external field, taken to 
be in the z direction, at some point the system will go 
unstable and have an avalanche. This will involve the 
nonlinear and possibly chaotic motion of its spins that 
interact through ferromagnetic and dipolar interactions. 

We will first analyze what happens during the 
avalanche when we set 7 to zero. This will describe 
the dynamics of the system for short time scales. In 
this case, the dynamics conserve energy. If the dynam- 
ics are sufficiently nonlinear, we might expect the system 
to be ergodic, which means that its equilibrium behav- 
ior is well described by the microcanonical ensemble. In 
this case, the avalanched state is not static but is one of 
a system at a finite temperature, and we will see that 
this temperature can be quite high, even above the ferro- 
magnetic transition temperature. However there are two 
reasons why ergodicity may break down. The first being 
if the Hamiltonian has strict rotational symmetry about 
an axis leading to conservation of angular momentum [7| 
in that direction Q. 

The second reason why ergodicity may be broken is 
more subtle and to our knowledge, this is the first time 
a scenario of this kind has been proposed. After the 
system initiates an avalanche, energy is transmitted into 
neighboring spins, some of which will be in the form of 
spin waves. This will propagate energy away from the 
avalanche region which will decrease the temperature of 
the avalanched spins, implying that as time progresses, 
the avalanched region becomes cooler. So if neighbor- 
ing spins are not recruited, the avalanche will be extin- 
guished. In this sense, the spin waves act as a damping 
term even if there is no damping in the LLG equation. 



For long times, still assuming 7 = 0, the energy will be 
distributed in all the degrees of freedom of the entire sys- 
tem, which means in the limit of infinite system size, the 
temperature of the system will have dropped back down 
to zero. One is then left with a system that has produced 
only a sub-system-size avalanche and has got trapped in 
another local minimum. 

Therefore we propose that an avalanche that is initially 
large enough, will propagate through the whole system 
causing it to go into a state of statistical mechanical equi- 
librium, often at a high temperature (for 7 = 0). If 
the initial avalanche is small, the avalanche will usually 
die out instead, leading to only a finite number of spins 
changing the sign of Sz ■ 

We now turn to two dimensional numerical experi- 
ments to support these claims and study the case of finite 
damping. Most real experiments on avalanche dynam- 
ics have been effectively two dimensional j^. Dipolar 
forces were not included as they complicate the analysis 
by adding an additional parameter. Their effects will be 
the subject of future work. 

We consider a Hamiltonian that couples nearest neigh- 
bor spins on a two dimensional square lattice and con- 
tains an anisotropy term where the orientation of the easy 
axis is randomized slightly about the z-axis and there is 
disorder in the ferromagnetic coupling. 
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where the J^'s are the ferromagnetic coupling constants 
drawn from a uniform distribution with nonzero positive 
mean, a is a measure of the anisotropy. We choose the 
fii to be random but biased towards the z-axis (out of 
plane). See ref. [3 for details. This system is placed in 
an external field B^xt- 

The system was started at high field and the field was 
lowered adiabatically by evolving the system at fixed Bext 
until, to a high accuracy, there was no further change in 
spin variables. To obtain convergence the damping was 
made finite, 7 = 1. After this the field was lowered again. 
An avalanche was defined to occur when the maximum Sz 
among all the spins changed by a finite amount A = 1. 
At that point, a successive approximation scheme was 
initiated to find the precise field at which the transition 
takes place, further evolution can then proceed using the 
same procedure. When an avalanche of desired size was 
detected, the system was restarted with the same exter- 
nal field and pre-avalanche configuration but now with a 
different value of damping and the evolution of the sys- 
tem was recorded. 

A common scenario is to find that the whole system 
will avalanche for sufficiently low damping, but will have 
a sub-system-size avalanche when the damping is above 
some critical threshold, that depends on the precise con- 
figuration right before the avalanche. This is demon- 
strated for a 128 x 128 spin system in Fig. Q] At 7 = 0.9, 
the post-avalanche cluster of avalanched spins has an ap- 
proximate diameter of 8 lattice spacings. However at 




FIG. 1: Snapshots of systems with different damping 7 start- 
ing with the same configuration, during an avalanche, (a) For 
7 = .8 and (b) For 7 = 0.01. The intensity represents the local 
spin motion. 



7 — 0.8, the whole system avalanches. Fig. ^ shows 
greyscale images of the motion of the system for 7 — 0.8 
(a), and 7 = 0.01 (b), during their avalanches. The inten- 
sity is proportional to \ds/dt\. In both cases, the system 
was started in the same configuration right before the 
avalanche. The motion of the system is confined to the 
cluster's surface for 7 = 0.8 but is spread out for 7 = 0.01 
in a ring-like structure, in which the magnetic system 
has an elevated effective temperature. In the case of no 
damping, many islands in front of the avalanche's main 
boundary form and this elevated temperature range en- 
compasses the entire avalanche region. After completion 
of the system-size avalanche, the entire system continues 
to move indefinitely in an ergodic phase. 

From these observations, one expects the hysteresis 
loop to change as one varies the damping in these simu- 
lations. This has been verified directly. 

For 7 = 0, the highly chaotic phase is energy conserv- 
ing and if there are no other invariants of motion, then 
the system should be well described by a microcanonical 
ensemble. For a large system, this ensemble is well known 
to be equivalent to the canonical ensemble (at finite tem- 
perature) . We verified by direct simulation of 32 x 32 size 
systems that this was the case by measuring static and 
dynamic correlation functions, for example, (si{0) -Sjit)). 

The finite temperature of the avalanche is often higher 
than the transition temperature. It is high because the 
pre-avalanche metastable configuration has an energy 
much higher than the minimum T — state. When 
the system gets out of its trapped static configuration it 
therefore has a lot of excess energy. By considering the 
Ising model with low disorder and estimating the crit- 
ical external field for avalanches to take place, we find 
that the system can have an energy close to zero, which 
is consistent with our numerical results of a high post- 
avalanche temperature [lOj 

We next turn to the mechanism by which this ergodic 
region spreads. As mentioned above, the effective tem- 
perature of the ergodic phase is quite high, with large 
amplitude motion over a short time-scale. The pre- 
avalanche configuration is static and when these two re- 



gions are connected together, there will be energy trans- 
fer between them. One would expect that over a large 
scale, Fourier's law should hold, so that that the temper- 
ature in the ergodic region will heat up the metastable 
region. Given thermal energy, the metastable region now 
has the opportunity to tunnel into the stable phase. 

We will now construct a simple one dimensional model 
that attempts to capture the above physics. We use a 
variable (j)i to denote if site i is part of an ergodic region, 
4>i = 1, 01 metastable region, (pi = 0. There is a temper- 
ature field, which starts off being zero in the metastable 
region and a non-zero constant To for the ergodic sites 
that seed the avalanche. This temperature corresponds 
to the energy released per spin when it becomes part of 
the ergodic region. The equation 



dT 
'di 



= DV^T + Tf 



dt 



-vT 



(3) 



describes thermal diffusion with diffusion coefficient £>, 
but adds a source term when the region becomes er- 
godic. In this case, V^ is a discretized second deriva- 
tive, oc Ti^i — 2Tj -|- Ti_i. We have also included for the 
sake of generality, a last term, vT, which is related to the 
damping in the system. This gives a time scale for the 
temperature to die out. For example in fig. QJa), large 
V corresponds to high temperature only on the surface 
of the cluster, whereas for low i^,^b), it persists over a 
fairly thick surface layer. However the case of i/ = (i.e. 
no damping) will be the focus of study below. 

The probability per unit time that ipi will go from 
to 1 is r(i), which we can take for example to be of the 
Arrhenius form ta = i^oexp(— l/T^), where Ti is the (di- 
mensionless) temperature on site i. 

Together these define a simple model for the disappear- 
ance of the metastable region. We are now in a position 
to analyze under what circumstances an avalanche prop- 
agates and when it dies out. We find that for sufficiently 
large Tq and initial width wq, of the ergodic seeded re- 
gion, the avalanche propagates indefinitely, but dies out 
if these two quantities are too small. 

As the avalanche propagates, the temperature at the 
surface will be Tq implying that the temperature in the 
interior will be almost constant with the same value. As 
the temperature diffuses to sites with = 0, there is a 
finite probability of them tunneling to (j) = \. 

To estimate the boundary as a function Tq and width 
w, we assume that for small T that r{T) is a rapidly 
increasing function of T. If the initial ergodic region 
is a top hat function at temperature Tq, then for large 
w, there is a long time when the temperature field next 
to the boundary will approach Tq/2. This time r will 



scale as 



Therefore the probability p that a site 



next to the boundary will change to </) = 1 should equal 
r(To/2)w^ [13I for p « 1. If it does not succeed in 
tunneling during this time r, the avalanche will die out, 
otherwise it will continue to propagate. As the avalanche 
spreads, w increases meaning that p increases so that it 
is easier to seed new sites. Therefore we expect the re- 
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FIG. 2: (Color Online) The boundary between avalanche 
propagation and extinction for the one dimensional model 
discussed in the text. The "-I-" symbols are the numerical 
values and the solid line is an analytical fit to the data. 



quirement for an avalanche to propagate is that 
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ri/2(ro/2) 
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where c is a constant. We have checked this numerically 
for the case of Arrhenius tunneling function r = ta^ 
as mentioned above. The results are shown in fig. El 
The plus symbols are the 50 percent probabilities of an 
avalanche propagating indefinitely, below the line, they 
will die out. The solid line is a best fit for the constant 
c in Wo = cexp(l/ro)- Given the simplicity of the esti- 
mate, the agreement is excellent. 

Local one dimensional models of avalanches with ran- 
domness will often not show a transition to propagation 
because there is a finite probability that at some time it 
will encounter conditions causing extinction. The differ- 
ence here is that the temperature field widens with in- 
creasing time. So if a site adjacent to the ergodic region 
fails to tunnel, in the mean tunnel time, the tempera- 
ture field will take longer to die away with increasing w, 
and therefore the probability of extinction rapidly goes 
to zero as the width increases, again for v = Q. 

Therefore for a zero temperature magnetic system with 
no spin damping, we expect that an initial disturbance 
will likely propagate if its size is above some threshold 
value, causing a transition between a static configuration 
and one at some effective finite temperature. Of course, 
if there is any coupling to degrees of freedom other than 
spin, so that 7 > 0, the motion will die out and the spins 
will stop moving. In such a case, the avalanched spins 
will lose their energy and cool down to zero temperature. 
The system is locally being annealed at a finite rate. 

One might expect that the inclusion of more realis- 
tic non-relaxational dynamics would not influence critical 
properties, because the ergodic region has a flnitc length 
scale for 7 > 0. If 7 is small, one would expect these ef- 
fects to shrink the critical regime. However there is also 
a possibility that the critical behavior would be altered, 
and this deserves closer scrutiny. Work on spin systems, 
with a conserved order parameter (e.g. the Heisenberg 



model without disorder) has shown that precessional mo- 
tion is relevant to the dynamic critical behavior J^ in 
equilibrium. Also, it is interesting to note that the few ex- 
periments available on hysteresis criticality do not seem 
to agree quantitatively with Random Field Ising Model 
calculations based upon relaxational dynamics [l3|- 

We now consider the issue of RPM. For the proof of it 
to be valid, a no passing rule must be satisfied 2] related 
to earlier work on charge density waves [iJl ■ If there are 
two configurations with the spins in the first all more neg- 
ative than in the second, then continued evolution under 
the same field will not guarantee that the second set of 
spins will remain more negative, as is required by this no 
passing rule. This is because the motion of the spins is 
effectively thermal and fluctuations can flip a spin in the 
second system above the value of the first system. This 
is very different than the relaxational dynamics needed 
to give RPM, where this can never happen. On the other 
hand, over a large enough scale, it may be unlikely that a 
coarse grained variable will violate RPM but nevertheless 
it is possible for this to occur. 

For the simplified model described above eqn. 13 finite 
damping i/ > can be analyzed. In this case it is quite 
similar to heat balance models for explosive crystalliza- 
tion [i^, which show many interesting properties |lall3- 
Instead of the temperature field at the front widening in- 
definitely, it should be of finite width, leading to a finite 



probability, per unit time, of the avalanche dying out. 
Thus one expects that in one dimension, propagation will 
always terminate eventually. In higher dimensions, be- 
cause the surface area of the front is increasing, we do 
expect to see infinite sized avalanches. 

Their are many interesting further questions to inves- 
tigate. If the total perpendicular angular momentum is 
only weakly broken by interactions, then what effect does 
this have on the dynamics? How does precessional mo- 
tion effect the size of the critical region of avalanche dy- 
namics? What is the effect of including disorder on the 
simplified model, eqn. O* And last, can these consid- 
erations be extended to better understand avalanches in 
granular media? 

In conclusion, we have shown that avalanches and hys- 
teretic behavior in spin systems are strongly infiuenced 
by precession and the strength of damping, where small 
damping makes a system more prone to large avalanches. 
However even with no damping term, coupling to spin 
waves can lead to the termination of avalanches. For fi- 
nite but small damping, the growth front becomes spread 
out over a large region, for which the spins inside can be 
described, for short times, by an ergodic system at high 
temperature and, for longer times, slowly anneal to a low 
temperature state. 

J.M.D. thanks A. P. Young and Trieu Mai for useful 
discussions. 
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